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Abstract 

Neutrino reaction cross-sections, (z^,^ - ), (u e ,e~), /i-capture and photoab- 
sorption rates on 12 C are computed within a large-basis shell- model frame- 
work, which included excitations up to 4hL0. When ground-state correlations 
are included with an open p-shell the predictions of the calculations are in 
reasonable agreement with most of the experimental results for these reac- 
tions. Woods-Saxon radial wave functions are used, with their asymptotic 
forms matched to the experimental separation energies for bound states, and 
matched to a binding energy of 0.01 MeV for unbound states. For compari- 
son purposes, some results are given for harmonic oscillator radial functions. 
Closest agreement between theory and experiment is achieved with unre- 
stricted shell-model configurations and Woods-Saxon radial functions. We 
obtain for the neutrino-absorption inclusive cross sections: a = 13.8 x 10 
cm 2 for the {y^,fi~) decay-in-flight flux in agreement with the LSND datum 
of (12.4 ± 1.8) x 10~ 40 cm 2 ; and a = 12.5 x 10~ 42 cm 2 for the (v e , e") decay- 
at-rest flux, less than the experimental result of (14.4 ± 1.2) x 10 -42 cm 2 . 
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I. INTRODUCTION 



One of the key neutrino-nucleus reactions measured at the two accelerator based neutrino- 
oscillations searches, LSND at the LANSCE facility in Los Alamos and KARMEN at the 
ISIS facility in the U.K., is the scattering of neutrinos from carbon in the liquid scintillation 
detector. The neutrino source at both experiments comes from the decay of pions produced 
in the beam stop. The vast majority of pions decay at rest, and the electron neutrinos 
thus produced have enough energy to cause a nuclear charge-exchange reaction on carbon, 
12 C(z/ e , e~) 12 N. At LSND 3.4% of the pions decay in flight, producing muon neutrinos of 
sufficient energy to interact via the reaction 12 C(z/ At , /x~) 12 N. The signal for the oscillation of 
decay-in-flight (DIF) muon neutrinos {y^ — > v e ) at LSND is the appearance of high energy 
electrons from the v e C — > Ne~ reaction. Extracting oscillation parameters from this search 
requires knowledge of the expected cross-section. In addition, the measured z^C^N/x - 
cross-section acts as a test of the DIF flux and of the detector efficiency The KARMEN 
experiment also has the efficiency to measure the vC — > Bfi + reaction. At both LSND 
and KARMEN the inclusive u e C— >Xe~ and the exclusive u e C—>- 12 N g s e~ cross-sections are 
measured for the electron neutrinos from the decay of the pion at rest. The flux of the 
decay-at-rest (DAR) electron neutrinos is given by the Michel spectrum, and the u e C cross 
sections provide a strong constraint on the nuclear-structure models used for carbon. 

The inclusive 12 C(z/ At , /i~)X cross section for the DIF flux at LSND has been measured. 
The first calculations for this cross section were carried out by the the Caltech group |TJ 
using a continuum Random Phase Approximation (RPA) model. The calculated cross sec- 
tions over-estimated experiment by almost a factor of two, and suggested that the measured 
cross sections may be inconsistent with other observables for 12 C. The most recent version of 
this work || allows for a partial occupancy of the nuclear subshells in the continuum RPA, 
a feature that brings theory and experiment closer together, but still does not remove the 
discrepancy fully. The other key observables that need to be considered are /x-capture, the 
DAR 12 C(z/ e , e~)X, (e, e'), photo-absorption, and /5-decay. These different probes involve 
different energy and momentum transfers and, thus, constrain different aspects of the cal- 
culations. The disagreement between theory and experiment for the DIF v^C cross section 
prompted new calculations ||(|] in an effort to uncover possible shortcomings of the RPA 
calculations. In the present work we perform a series of shell-model based calculations, which 
include excitations up to 4hu) for 12 C, 12 N and 12 B, each involving different assumptions and 
appr oximat ions . 

The DAR neutrino spectrum involves neutrino energies — 52 MeV, with an average 
neutrino energy E v ~ 32 MeV. The Q- value for the charge-exchange reaction (z/ e , e~) on 12 C 
is about 17 MeV. The DAR inclusive cross-section is then dominated by low multipoles (1 + , 
l - , 2~) and by excitation of the giant dipole resonances. In the case of the u^C cross section, 
the DIF muon neutrino flux involves an average neutrino energy of about E v ~ 150 MeV, 
but the flux is finite up to E v ~ 250 MeV. The Q-value for the (V M , fi~) reaction is close to 
123 MeV. Calculations for this cross section need to include both a good description of the 
giant resonance region (E x ~ 15 — 40 MeV of excitation in 12 C) and of higher excitation 
energy regions (up to 80 — 100 MeV). Furthermore, all multipoles A ~ — 5 make significant 
contributions to the inclusive cross section. 
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II. GENERAL NUCLEAR STRUCTURE CONSIDERATIONS 



In the simplest shell model, 12 C consists of four neutrons and four protons in the p- 
shell outside a closed 4 He core. Excited states reached by the (u, £~) reactions are simple 
particle-hole states built on the ground state. However, the structure of both the 12 C g . s . and 
the continuum states of 12 N also involve more sophisticated configurations and configuration 
mixing. There is a limit on the size of the model space that can be included in any calculation, 
so that calculated cross sections necessarily involve some level of approximation. In the case 
of the neutrino reactions on carbon the approximations that are most likely to affect the 
predicted cross sections are (a) the 12 C ground state p-shell structure, (6) the treatment of 
ground state correlations beyond the p-shell, (c) the model space truncation, especially for 
the final states, (d) configuration mixing in the final states, and (e) the nuclear radial wave 
functions. In this paper we examine the effect of each of these on the predicted cross section. 

An approximation that is inherent in the continuum RPA calculations of Kolbe et al. 
|]J is the restriction of the 12 C ground state to a closed 0p 3 / 2 shell. This approximation is 
necessary to build a spectrum of particle-hole states representing the excitations in 12 N and 
to evaluate their transition cross sections in the RPA. That this approximation is very poor 
for the lowest positive-parity states of configuration (p^^Pi/2) is we U known and recognized 
in the calculations of Kolbe et al. 0, who reduce their calculated cross section to the 12 N ff S . 
by a factor of four. A smaller reduction factor of order 1.5 is obtained when this work is 
extended [|| to include partial occupancies for the pi/2 subshell. An additional reduction 
of the cross section to the other low-lying positive parity states, particularly to the 2 + 
state at 0.96 MeV of excitation in 12 N, should also be included. It is unclear whether any 
further correction is required for particle-hole states of energy lhu, 2Hu . . . above these 
lowest positive-parity states in 12 N. Kolbe et al. have argued against additional suppression 
factors, other than the one that has been applied to the GT transition to the 12 N ground 
state. 

In the case of particle-hole excitations out of the p-shell, the p-shell structure of the 
ground state does not affect the total sum-rule for a given operator, which is determined 
by the number of particles in the p-shell. The total spin independent AS = multipole 
strength is independent of the p-shell structure of the ground state. However, the total spin- 
dependent AS = 1 strength is distributed differently over the different spin-multipoles for 
different p-shell ground states ||^| . Thus, for flux-averaged neutrino reactions the predicted 
cross section will, in general, be model dependent, reflecting the model ground state spin 
structure. 

A reasonable approximation for the structure of 12 C is the p-shell equivalent of three 
a-particles. This corresponds to an L = S = ground state, with good SU(4) symmetry, 
[444]. The Cohen-Kurath interaction ]7) predicts that this state makes up 78% of the 12 C 9 . S . 
wave function. In contrast, the closed p 3 / 2 -state contains 16% S = and 6%[444] symmetry, 
Table |. Thus, the spin response for the closed P3/2 state will differ from that of the Cohen- 
Kurath ground state. Determining the degree to which this model dependence is reflected 
in the predicted neutrino cross sections is a main aim of the present work. 

Proton threshold in 12 N is at 0.601 MeV, so that all levels but the ground state are 
unstable. States up to about 8 MeV of excitation in 12 N have been studied via the 12 C(p, n) 
and (n, p) charge exchange reactions. The structure of these states can be understood largely 



3 



within a (0+l)hu shell-model calculation, although the predictions of these Tamm-Dancoff 
calculations overestimate the charge exchange cross sections by about a factor of two. Little 
detailed information exists for higher energy regions of the 12 N continuum. 

The unbound nature of the states in 12 N is in strong contrast with the deeply bound 
p-shell particles in the 12 C ground state. Thus, there will be a strong lack of overlap between 
the initial- and final-state radial functions. The continuum RPA calculations of Kolbe et al. 
|]J explicitly treat the continuum nature of the excited states of 12 N. In contrast, the present 
shell model calculations are discrete state calculations. Calculations of the 12 C(p,n) and 
12 C(n, p) reactions which took binding energy effects into account via the use of Woods-Saxon 
single-particle wave functions found a ~30% suppression of the predicted cross sections over 
the predictions using harmonic oscillator wave functions. In the present work we believe it 
is essential to work with Woods-Saxon radial functions, but for comparison purposes we will 
give some results with harmonic oscillator single-particle functions. 



III. SHELL-MODEL CALCULATIONS 
A. Model Space 

To investigate the effect of the approximations inherent in any model calculations of the 
inclusive neutrino cross sections on carbon we mount a series of shell-model calculations, 
and present calculations in four separate model spaces. The spaces include excitations up 
to 4hu, and we label these model spaces (1) Closed-shell TDA, (2) Closed-shell RPA, (3) 
Closed-shell RPA +2p-2h, and (4) Unrestricted shell model. The labels represent the model 
spaces: 

1. Closed-shell TDA. The 12 C ground state is described as a closed 0p 3 / 2 shell and the 
12 N excitations as one particle-one hole (lp-lh) states: 

| 12 C> = |0>, 

\ 12 N) = \(h- 1 ,p) nn "J,T) n = 0,1, 2, 3, 4 (1) 

The lp-lh states are labeled by their energy of excitation, nh~u, in an oscillator model. For 
example, states ip^/ 2 ,Pi/2) are Ohu excitations, (s^} 2 ,Pi/ 2 ) an d (p^ 2 , sd) are ITiui excitations, 
(s^ 2 ,sd) and (p^ 2 ,pf) are < ^ IUJ excitations, and so on. We present calculations up to Atiui 
excitation. A shell-model calculation in this model space is equivalent to the Tamm-Dancoff 
Approximation (TDA) for particle-hole excitations. 

2. Closed-shell RPA. In this case, 2p-2h excitations are included in the 12 C ground-state 
wave function 

\ 12 C) = I 0)+ I {hi\h?)J x ,T x - ( Pl ,p 2 )J 1 ,T 1 : 00) 

\ l2 N) = \{h~\p) nhw J 1 T) n = 0,1, 2, 3, 4 (2) 

Here h±,h 2 span the hole orbitals, Osi/2 and 0^3/2, while pi,p 2 span particle orbitals 0p±/ 2 , 
0d 5 / 2 , lsx/2, . . .. The highest-energy orbital included in the particle space matches the 
highest-energy orbital, p, in the lp-lh basis of energy nhu. Although we have labeled this 
calculation, RPA, there are more 2p-2h states included in the 12 C wave function here than are 
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normally present in an RPA calculation. This is because the 2p-2h correlations introduced 
in RPA are restricted to the type 

| (hj\ Pl )^J,T; (h 2 \p 2 ) nh "J,T : 00), 

where the 2p-2h states are made up only from the coupling of two lp-lh states of spin, isospin 
J, T that comprise the basis states of 12 N. Not only is this 2p-2h basis smaller than that in 
Eq. (§) it also is not fully antisymmetrized. Hole states h± and h 2 , and particle states p\ 
and p 2 are not antisymmetrized with respect to each other. This shortcoming is not present 
in the basis of Eq. (0), where the coupling order shown makes it easy to antisymmetrize the 
states. Although, as we have just explained, a calculation in this model space is more than 
just RPA, this case is the closest we have to the calculations of Kolbe et al. [||. 

3. Closed-shell RPA + 2p-2h. In this case, we add some 2p-2h configurations to the 
lp-lh basis states of 12 N: 

\ l2 C) = | 0)+ | (h^ 1 ,h 2 1 )J 1 ,T l] ( Pl ,p 2 )Ji,T 1 : 00) 

| 12 iV) = | {h-\p) n ^J,T)+ | ( P 3 / 1 2 ,/i- 1 )^,T3;(p 1/2 ,p)J4,T 4 : J,T) n = 0,1,2,3,4 (3) 

Note that among the 2p-2h configurations for 12 N, one of the holes is restricted to the 0p3/ 2 
orbital and one of the particles to the 0px/ 2 orbital. In this way, we are moving further away 
from the assumption that the 12 C wave function is a closed ps/ 2 shell. Again the other hole 
orbital, h, may span the hole orbitals, 0si/ 2 and 0p 3 / 2 , while the other particle orbit, p, spans 
orbitals 0pi/ 2 , 0d§/ 2 , lsi/ 2 , ... up to the same maximum characterizing the lp-lh states. 

4. Fully Unrestricted Shell Model. Finally, we move to an unrestricted shell-model 
calculation, with full mixing between all configurations: 

\ 12 C g.s. > = I (P3/2P1/2) 8 : 00 > + I (p 3 /2Pi/2) 6 (ls0d) 2 : 00 > 

+ I (P3/2Pi/ 2 ) 7 ( W) : 00 > + I (0 S )- 1 (p 3 /2Pi/ 2 ) 8 (ls0rf) : 00 > 

| 12 AT 7T+ > = I (p 3/2Pl/2 ) 8 : JT > + I (p 3/2 p 1/2 ) 6 (ls0d) 2 : JT > 

+ I (P3/2Pi/2) 7 ( W) : JT > + I (0 S )- 1 (p 3 /2Pi/ 2 ) 8 (ls0d) : JT > 

\ 12 N 7T- > = I ( P3/2 p 1/2 ) 7 (ls0d) : JT > + I (Us)-\p V2 p 1/2 f :JT>. (4) 

Spurious center-of-mass states were eliminated exactly by making a transformation to an 
orthonormal basis for which the expectation value of the centre-of-mass Hamiltonian is zero, 
< H cm >= 0. Our 2%uj shell model basis contains ~2000 states for each multipole, which 
is reduced by ~ 25 when transformed to a non-spurious basis. Because of computational 
limitations we are only able to take our "unrestricted" shell model calculations up to 2Tnu. 
However, we assume an extrapolation of a 2hu shell-model calculation to 4hu scales in the 
same way as the closed-shell RPA + 2p-2h calculation. 



B. Effective Interactions 

Having selected the appropriate model spaces, we must specify the operative effective 
interactions to use in these spaces. We choose interactions that have been fitted to repro- 
duce experimental spectra, namely: For interactions among the p-shell orbitals, we use the 
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interaction of Cohen and Kurath (8-16)2BME; for interactions among the s, <i-orbitals, 
we use the universal s, d- interaction (USD) of Wildenthal ||; for interactions among the 



p, /-shell orbitals, we use the G-matrix of Kuo and Brown || as modified by Zuker |L0| and 



known as KB3. All other two-body matrix elements, including the cross-shell interactions, 



were calculated using the Millener-Kurath [11] parameterization in terms of central, spin- 



orbit, and tensor interactions, whose radial forms are Yukawa functions of various strengths 
and ranges in different spin-isospin channels. The parameters are determined in fits to spec- 
tra of unnatural parity states in p-shell nuclei, and so are particularly appropriate for the 
cross-shell interactions between the p- and s, <i-shell shells. We, however, use this interaction 
for all cross-shell matrix elements required. 

Finally, the single-particle energies must be specified. For the (8-16)2BME, USD and 
KB3 effective interactions these energies are given as part of the fitted interaction, but are 
referenced to their respective cores in their shell-model usage, namely 4 He, 16 and 40 Ca. 
In the present work, we wish to make our reference core a closed p%/2 shell at 12 C, so these 
single-particle energies are shifted according to the formula 

4 = 4 + E {2J to ] % +1) ^ h ^ T I v I ^) J ' T >' ^ 

h,J,T Z \ Z J L ) 

where ef is the single-particle energy of an orbital, j, relative to a core A, ef relative to 
a core B, A > B, and the sum, h, is over all the states occupied in A but unoccupied in 
B. Since only the relative separation of single-particle energies is relevant in a shell-model 
calculation we set the single-particle energy of the Opi/2 orbital at its experimental value 
of —4.95 MeV. For orbitals above the p, /-shell, we have little guidance in choosing the 



single-particle energies. We use, therefore, a formula given in Bohr-Mottelson [12 



const. + vi s hw(l.s) + vuhu (l 2 - (1 2 )at) , (6) 



(\ 2 ) N = 1 -N(N + 3) 1 (7) 

where N = 2n + / is the principal quantum number for the oscillator orbital, n the number 
of radial nodes (excluding the origin and infinity) and I the orbital angular momentum 



quantum number. We use the values vi s = —0.127 and vu = —0.03 from [ji^j . The choice of 
Tiuj to be used in this formula should be one appropriate for nuclei where the (s, d, g)- and 
(p, /, /i)-orbitals are the valence orbitals, not an huo appropriate for 12 C. We choose tku = 7.2 
MeV. The constant in Eq. (||) is again chosen so that when the formula is used for the Opi/2 
orbital it reproduces the experimental value, —4.95 MeV. The values of the single-particle 
energies used are given in Table |I|. 

IV. FORMALISM 

We apply the results of shell-model calculations just described to the evaluation of a 
number of isovector weak and electromagnetic transitions in the A = 12 system. We consider 
the inclusive neutrino reactions 
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u^+^C^X + fT, (8) 

u e + 12 C -> X + e~ (9) 
where X is unobserved, the muon capture process 

M -(15')+ 12 C^X + ^, (10) 

and photoabsorption 

7+ 12 C^X, (11) 

which proceeds only through the isovector El multipole. In addition, we take guidance from 
earlier calculations [p~3| , |T4|| for the (p,n), (n,p) and (e, e') reactions to the isovector giant 



resonances of interest. Further we consider the exclusive cross section to the ground state of 
12 N for the first three of the above reactions, as well as the Gamow- Teller beta decay from 
12 N to the ground state of 12 C. 

The formalism for calculating these cross sections has been given by O'Connell, Donnelly 
and Walecka JT5fl . For the neutrino absorption reactions, Eqs. (|8|) and (^), the cross section 
is given by 

a(E p ) = ^Y. [ +1 d(cos6)5(E l -E f + E u ~E e )p e E E |M| 2 , (12) 

where G is the weak interaction coupling constant, G/(hc) 3 = 1.16639 x 1CT 5 GeV -2 , p u , E v 
are the neutrino momentum and energy in the lab system, p£, En are the outgoing lepton 
momentum and energy, cos# is the cosine of the angle between the electron and neutrino 
directions, and Ei — Ef is the mass difference between the initial and final nuclei. It is 
convenient to define an energy and momentum transfer: q = Ei — Ef = Ei — E v and 
q = pi — p u . There is a sum over all final nuclear states, which in our calculations becomes 
a sum over all the states available in a given shell-model calculation. 
Lastly \M\ 2 is given schematically 

oo 

I M | 2 = £(LiWi (A) + L 2 wi X) + L 3 W^ X) + L 4 wi X) + L 5 I^ 5 (A) ) (13) 

A=0 



where L± . . . L$ are the five lepton traces given in Table II of Ref. |y| and are functions of 



go, q, and cos 6, while W± X ^ . . . are certain combinations of squares of nuclear matrix 
elements. The sum is over all multipolarities of operators, A, satisfying angular momentum 
condition A(Jj, Jj, A), where Jj and Jf are the initial and final spins. 

In general, there are seven transition operators, 0^ x \qx.), detailed in the tables of Don- 



nelly and Haxton ||16|| . The operators are functions of q, the magnitude of the momentum 
transfer, q —\ q|, and x, the position coordinate. Reduced matrix elements of these operators 
can be decomposed into two factors: 

(J f ||0 (A) (gx) || Ji) = £ (J f || [al,a jf3 ]^ || JOOa H0 (A) O?x) || fr), (14) 
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where the accompanying isospin quantum numbers have been suppressed for economy of 
notation. The first factor is the expectation value of shell-model creation and annihilation 
operators evaluated with the many-body shell-model wave functions and is known as the one- 
body density matrix element (obdme). The second factor, the single-particle matrix element 
(spme), is independent of the many-body wave functions but depends on the transition 
process in question through the momentum transfer, q. This factorization enables the shell- 
model calculations to be separated from the transition rate calculations. Files of obdme 
prepared in the former are transferred to the latter. The spme also involve radial integrals of 
spherical Bessel functions, Jl{qx), or their derivatives with two single-particle nuclear radial 
functions. For the latter, we use Woods-Saxon functions. The nuclear matrix elements are 
also multiplied by form factors, which are functions of q, as given in Ref. [15| . 



In the shell-model calculations, the center-of-mass is treated on an equal footing with 
all other 3 (A — 1) coordinates, and which causes an erroneous center-of-mass contribution 
to the neutrino cross sections and muon capture rate. To correct for this we multiply the 
computed cross section by the Tassie-Barker function | <7 cm (g) | 2 , where 

9UQ) = e y/A (15) 

with y = (bq/2) 2 , q the momentum transfer, b the oscillator length parameter and A the 
nuclear mass number. For y ~ 1 and A = 12, this correction provides a 20% increase in the 
computed cross section.^] 

In the expression for the neutrino absorption cross-section, Eq. (|12"D, one further correc- 
tion has to be applied that represents the distortion of the outgoing lepton in the Coulomb 
field of the daughter nucleus. We are guided by the work of Engel |17| in this matter. The 
energy of the emerging lepton, while under the influence of a constant electrostatic potential 
within the nucleus, is effectively shifted: 



E £)CS = Eg — V(0) 



P£,es = [El cS -m 2 



1/2 



™=-if (i6) 

where Eg, pe and m are the energy, momentum and mass of the emerging lepton, Z the 
charge number for the daughter nucleus (with a plus sign for particle leptons, and a minus 
sign for antiparticle leptons), R the nuclear radius and a the fine-structure constant. If 
PiR is less than ~ 0.5, then the lepton wave function is primarily s-state, and the usual 
Fermi function, F(Z,Ei) is used as a multiplicative correction factor in Eq. (|12|). For DAR 
cross sections, piR is always less than 0.5. For DIF cross sections, the emerging muons 



have piR > 0.5 for most of the energy range. Then we follow Engel's suggestion |T7j that 
the cross-section be multiplied by a factor pi tCS Ei jeS /(p e Ei) and the momentum transfer 
q = |q| used in the spherical Bessel functions in the radial integrals be shortened by using 



1 This correction has been derived with the use of oscillator functions. We apply the same correc- 
tion even when using Woods-Saxon radial functions. 
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q e ff = pe,cs — Pi/- The multiplicative factor increases the DIF cross section by about 16% 
while the shortening of q in the Bessel functions reduces this enhancement to about 11%. 

In the results to be given in the next section, for the neutrino cross sections Eqs. (H) and 
(|), we will give the flux-averaged cross-sections defined as 



a = J dE v k{E v )o{E v )l J dE u k(E v ) 



(17) 



where A(E V ) is the incident neutrino flux. For the liquid scintillator neutrino detector 
(LSND) experiment |Tj|], we use their most recently determined muon-neutrino flux [TJ5| 
from pion decays in flight. For the experiments by the KARMEN Collaboration [20|, the 



electron-neutrino beam comes from muon decays at rest and the flux is given by the Michel 
spectrum. 

The muon capture rate, Eq. (|TU|), is given by 



A, 



G 2 El 
2tt 



Ml 



(18) 



where the outgoing neutrino energy is E v = nig — €b + Ei — Ef, with mg the muon reduced 
mass, €b the muon K-shell binding energy and Ei — Ef = qo is the mass difference between 
the initial and final nuclear states. The factor (1 + E v jmT)~ x is a recoil correction, with 
mx the target nucleus mass. Further, <f>\$ is the muon K-shell bound state wave function 
evaluated at the nucleus 



'15 lav = 



-{Zam.f 



(19) 



with R a reduction factor for the finite size of the nuclear charge distribution, R = 0.86 for 
12 C. The total inclusive capture rate is the sum of Eq. (0), summed over all final states 
contained in the shell-model calculation. 

Finally, the photoabsorption cross section, Eq. flTID, is given by |T2 



a(E^ 



167T 

~9~ 



aEr r B(El; i -> f)8{E 1 - E f + Ei 



(20) 



where E 1 is the photon energy and B(E1; i —> f) the reduced transition probability for El 
photon absorption. Again the total absorption cross section, <r tot , is given by the sum of Eq. 
(p0|) , summed over all final states contained in the shell- model calculation. In the different 
calculations, the total absorption cross section reaches a saturation value at a photon energy 
Ej of the order of 50 MeV, slightly less in ltkv calculations, slightly more in (l+3)Hu. Note, 
only the El multipole is retained in the calculations. The experimental total absorption 
cross section ||27||, however, shows no sign of saturating and keeps growing with increasing 



E 7 . This is because higher multipoles are contributing. In quoting an experimental value, 
we have cut off the contributions at E 1 = 50 ±5 MeV for a total cross section of er tot = 21 ±1 
fm 2 . 



V. RESULTS 
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A. Inclusive Cross-sections 



We present the results of our calculations in Tables fT| to |VT|. In most of our tables the 
tabulated results are for Woods-Saxon radial wave functions. In Table |V| we list the results 
for three sets of radial wave functions, Woods-Saxon, and harmonic oscillators with b = 1.82 
fm and b = 1.64 fm. For each of the four model spaces, we give the results as a function of 
the size of the particle-hole space, increasing from Ohu excitations up to 4hu. In Table U I 
are the inclusive cross sections for (z/^, /i~)-scattering, (u e , e~)-scattering, muon capture, and 
photoabsorption on 12 C. The calculation in the TDA approximation overestimates the cross 
sections in all cases, while the two variants of RPA reduce these cross sections by about a 
factor of two and are in closer agreement with experiment. For the cross section, 

the calculations to Ahu in the closed-shell RPA, closed-shell RPA + 2p-2h, and unrestricted 
shell model yield 18.2, 16.7 and 13.8 x 10~ 40 cm 2 , the latter value being an extrapolation 
from a 2tkv calculation. This latter value, also, is within range of the experimental value of 
(12.4 ± 1.8) x 10~ 40 cm 2 |I8" |. This is our principal result. In ref 0, with the p%/2 subshell 



taken as a closed shell, Kolbe et al. obtain a result 19.8(5) x 10~ 40 cm 2 , where the uncertainty 
represents a spread between different choices of the effective interaction. Recently [[|, this 
work has been extended to include a partial occupation of the pi/ 2 subshell obtaining a result 
some two units smaller of 17.6(2) x 10~ 40 cm 2 . The present calculation reduces this four 
units further with the use of fully unrestricted shell-model configurations and Woods-Saxon 
radial wave functions. 

However, in obtaining an improved result for the DIF cross section, we must check that 
the same calculation still satisfactorily reproduces other relevant data. In Table |H| we 
therefore give results for (v e , e~) DAR cross sections, //-capture rates and photoabsorption 
cross sections. For (i/ e , e~) neutrino absorption, we calculate in closed-shell RPA, closed-shell 
RPA + 2p-2h, and unrestricted shell-model cross sections of 21.9, 20.4 and 12.5 xlO -42 cm 2 
respectively compared with the KARMEN result of (14.0 ± 1.2) x 10~ 42 cm 2 and the LSND 
result of (14.8 ± 1.3) x 10 -42 cm 2 . Again the unrestricted shell-model calculation gives a big 
reduction to the neutrino absorption cross sections resulting, in this case, with a value that 
is on the low side compared to experiment. However, it is only outside the experimental 
range by one standard deviation. By comparison, the continuum RPA calculation of Kolbe 
et al. obtains a result 14.4 x 10~ 42 cm 2 exactly in the right range. Thus it would seem that 
in obtaining a smaller DIF cross section as required, we are simultaneously underpredicting 
the DAR cross section. 

For /i-capture and photoabsorption, our results in Table [TTI| indicate the unrestricted 
shell-model calculation slightly overpredicting the experimental result, but only by about 
10%. In summary, then, our improved result for the DIF cross section comes with some 
modest deterioration in the DAR, //-capture and photoabsorption results. 



B. The Radial Wave Functions 

In discrete-state shell-model calculations with effective interactions, the choice of radial 
wave functions to use in evaluating transition matrix elements remains unspecified. There is 
therefore some freedom in making this choice. In Table [IV] we give results for three choices 
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of radial functions for the inclusive cross sections, but in this case explicitly exclude the 
ground-state contribution. 

A description of the rms charge radius of 12 C using harmonic oscillator single-particle 
wave functions requires an oscillator parameter of b = 1.64 fm. However, to obtain rea- 
sonably good fits to the shape of the electron-scattering form factors in 12 C considerably 
different values of b are often needed. An analysis |H],|n|[l4]] of the (e, e') form factors and 
(p,p f ) inelastic scattering data suggests the need for state dependent oscillator parameters 
in the broad range 1.64 — 1.94 fm. The problem arises as a result of the inability to describe 
loosely bound or unbound states of 12 N using the same oscillator parameter as is needed to 
describe the deeply bound nucleons in the 12 C ground state. 

In a systematic study of the effect of more realistic single-particle wave function on 
the 12 C(p, n) reaction Millener et al. |T3|JT4]] concluded that the use of Woods-Saxon wave 



functions reduces the predicted cross sections by about 30% over the harmonic oscillator 
predictions. Ohnuma et al. drew a similar conclusion. The large suppression of the 
cross section with the use of Woods-Saxon single particle wave functions is because of the 
large mismatch between the radial wave functions describing the l2 C g , s , versus the 12 N states. 
An effect of similar magnitude was found for the 12 C(e, e') form factors. 

The cross sections obtained using Woods-Saxon wave functions have been found PTJJ2BI 
to be quite similar to those obtained when the harmonic oscillator parameter is chosen to 
fit best the (e, e') form factor. Brady et al. [13|] obtained a good fit to the (e, e') form factor 
for the 12 C S . S . — > 12 C(15.11 MeV) using an oscillator parameter b = 1.82 fm. To examine the 
sensitivity of the predicted neutrino cross sections to the assumed oscillator parameter we 
recalculated the cross sections using a smaller value of b, and the results are tabulated in 



Table |V[ Relative to the b = 1.82 fm calculation the inclusive (z/ M , /z") DIF cross section 
increased only slightly for b = 1.64 fm, while the (V e , e") DAR cross section decreased by 17% 
and the //-capture rate decreased by 14%. The larger value of b in all cases gave predictions 
closer in agreement with experiment. 

The difference in the response of the three reactions to the larger oscillator parameter 
reflects the different average momentum transfer involved in each, DIF (q ~ 1 fm -1 ), DAR 
(q ~ 0.2 fm -1 ), and //-capture (q ~ 0.5 fm -1 ). Changing the oscillator parameter changes 
the shape of the axial-vector and vector form factors entering the neutrino cross sections. 
For momentum transfers below the first maximum of the form factor the neutrino cross 
section increases as b increases. In contrast, momentum transfers beyond the first peak in 
the form factor result in the predicted cross section being decreased as b increases. 

We next consider the use of Woods-Saxon radial wave functions. The mismatch between 
the initial and final radial wave function, resulting from binding energy changes, cannot be 
accounted for within a harmonic oscillator basis. It is impossible to treat the binding energy 
of the final nucleon rigorously within a discrete shell model, but an estimate of the effect of 
binding energies on the predicted cross sections can be obtained using Woods-Saxon wave 
functions that are only just bound. The Woods-Saxon result reduces the predicted inclusive 
DIF (z/ M ,/i - ) cross section by about six units, putting the result within the experimental 
range. However the predicted DAR (z/ e ,e") cross-section is similarly reduced by 2.8 units 
and is too small compared with experiment. The calculated /i-capture rate increases and 
slightly overpredicts experiment. 

The use of Woods-Saxon functions causes some radial integrals to become non-zero, 
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where they are identically zero for oscillator wave functions. This is seen easily for the 
El 7-ray transition matrix elements. With oscillators, the predicted B(E1) strength in 
(0 + 1 + 2 + 3)hu is slightly less than in a (0 + l)Hu) calculation. But with the use of 
Woods-Saxon wave functions the (0 + 1 + 2 + 3)hu> calculation is about 40% larger. This 
is due to the non- vanishing integral / R(n, l)rR(n', l')r 2 dr, when radial functions, R(n, I) 
and R(n',l'), represent single-particle states differing in energy by 3tku. The same effect is 
noticeable in the calculated /i-capture rate and to some extent in neutrino scattering cross 
sections. 

C. The Exclusive Cross Section to the 12 N„ , 

Next, we examine in detail the exclusive cross section to the 12 N ground state. The 
results are in Table [V]. The first row of the TDA calculation labeled Ohu is the single- 
particle + —> (p^ 2 ,pi/ 2 )l + transition rate, while the first row of the unrestricted shell-model 
calculations represents the complete p-shell mixing as given by Cohen-Kurath (8-16)2BME 
wave functions. We note that the cross sections in the configuration-mixed calculations are 
a factor four to five smaller than the single-particle estimate. The closed-shell RPA variants 
are only able to obtain about a factor of two of these reductions. Increasing the particle- 
hole space to include 2hu> and 4hu excitations only has an impact in the 10 — 20% range for 
neutrino absorption, muon capture and beta decay. 

In our complete unrestricted (0 + 2)hu calculations we paid particular attention to han- 
dling the strong Ahu = 2 interaction in a consistent way. When empirical or realistic 
Hamiltonians which give a reasonable description of Ohoo and ltkv spectra are used in full 
(0+2)huj spaces, the Ahu = 2 interaction pushes the ground state down by several MeV. 
The resulting spectrum is in poor agreement with experiment. With no restrictions placed 
on our (0+1 +2) fro; calculations we found that the giant monopole resonance was predicted 
to lie too low in the spectrum and large percentages of lplh excitations were predicted in 
the ground state wave function. To correct for this pathological behavior we adjusted diag- 
onal matrix elements of the Hamiltonian so as to restore the energy of the giant monopole 
resonance to lie above 20 MeV, as expected by self-consistent RPA calculation using density- 
dependent interactions. In addition we set the strong (lplh 2hu\H\0huj) matrix elements, 
that transform as (A,/i) = (2,0) under SU(3), to zero. The resulting 12 C ground-state wave 
function yields a predicted strength for the GT transition to the 12 N S s that is in reasonable 
agreement with experiment, Table |V]. 

D. The Multipole Decomposition 

Finally, in Table [VT] we give the breakdown in the contribution to the cross sections 
from each multipole, for oscillator wave functions using b = 1.82 fm. The most striking 
result is that 80% of the (u e , e~) inclusive cross sections comes from the ground-state 1 + 
multipole, with the remaining 20% coming from the excited state 1 + , 1~, and 2~ multipoles. 
Muon capture also is dominated by the ground-state 1 + multipole, being, in this case, 65% 
of the inclusive rate. The remaining 35% is dominated by the 1~ multipole with lesser 
contributions from the 2~, 0~, excited 1 + , and 2 + multipoles. By contrast, the (z/^,/i~) 
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inclusive cross section gets only a 10% contribution from the ground-state 1 + multipole, 
while the remaining 90% is spread over many multipoles, 1~ to 4r and 1 + to 3 + dominating. 

VI. CONCLUSIONS 

We have examined the (i^,// - ), (v e ,e~), and /i-capture reactions on 12 C within a shell- 
model calculations that included up to Ahcu of excitation. Calculations were carried out 
for four separate model spaces, for three sets of radial wave functions. Closest agreement 
between theory and experiment was obtained for the most sophisticated model spaces, which 
included both RPA-type correlations and an open p-shell 12 C ground state. 

Previous analyses of the (e, e'), (p, n) and (n,p) reactions have shown that when the large 
difference in binding energies between the 12 C ground state and the 12 N final states in taken 
into account, the predicted cross sections are reduced by about 30% relative to harmonic 
oscillator shell-model results. To examine the effect of more realistic single-particle wave 
functions we compared the predictions using two oscillator parameters, b = 1.64 fm and 
b = 1.82 fm and Woods-Saxon wave functions. The lower value of b is suggested by the 
12 C ground-state charge radius, while the larger value comes from a fit to the shape of 
the transverse form factor for the 12 C(0 — > 15.11; 1 + ) transition. Increasing b increases the 
predicted cross section for all low momentum transfer reactions. The photoabsorption cross 
section increases with b 2 , and the DAR (u e ,e~) cross-section increased by a smaller but 
similar amount. For higher momentum transfer processes (q ~ 1 fm" 1 ), such as the DIF 
neutrino reaction and the (p, n) reaction at intermediate energies, the predicted cross section 
decreases with increasing b. Of course, binding energy differences result in a mismatch 
between the initial and final radial wave functions which cannot be mimicked by a change 
in the oscillator parameter, and fits to the (e, e') form factors for states in 12 C within an 
oscillator basis find a need for state-dependent oscillator parameters in the range 1.64— 1.94 
fm. Incorporating binding energy effects through the use of Woods-Saxon wave function has 
a large effect on the predicted cross sections for all reactions. The predicted inclusive DIF 
(i^,// - ) cross section is in agreement with the experimental range. However, the Woods- 
Saxon calculation underpredicts the DAR (z/ e , e~) cross section and overpredicts the //- 
capture rate. 

Given the model dependence noted in the present set of calculations, it is difficult to 
determine the expected DIF cross section to high accuracy using constraints from other 
observables. For example, the ratio of the predicted DAR to DIF cross sections changed 
by about 20% for the two different values of oscillator parameter considered. Furthermore, 
we were unable to include a fully realistic treatment of the continuum nature of states of 
12 N. However the success of the Woods-Saxon calculation for the DIF cross section while 
only failing to obtain the experimental result for the DAR cross section, /i-capture rate and 
photoabsorption cross section by one to two standard deviations indicates the importance 
of a continuum treatment for the final states, coupled with an open-shell description of the 
nuclear configurations. 
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TABLES 



TABLE I. Spin structure of the 12 C ground state 



Model 








S = 




S = 1 




S = 2 




%[444] 


Three a- 


-particles 






100% 














100% 


Cohen-Kurath Interaction 




81% 




18% 




<2% 




78% 


(P3/2) 8 








16% 




59% 




25% 




6% 








TABLE II. Sinj 


de-particle energies. 










Osi/2 


0f»3/2 


Opi/2 


04/2 


1*1/2 


04/2 


O/7/2 


1P3/2 


O/5/2 






-34.40 


-12.26 


-4.95 


5.98 


3.70 


8.97 


10.16 


8.78 


9.39 


8.39 






14/2 


2si/ 2 


Oc/7/2 


14/2 


O/in/2 


I/7/2 


2p 3 /2 


Ohg/2 


I/5/2 


2^1/2 


13.87 


18.09 


20.46 


18.28 


20.54 


20.17 


25.32 


28.61 


25.56 


28.75 


30.08 
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TABLE III. Inclusive cross sections involving a continuum of 12 N states, including the ground 
state. Woods-Saxon radial functions are used, with their asymptotic forms matched to the exper- 
imental separation energies for bound states, and matched to a binding energy of 0.01 MeV for 
unbound states. 

(ffi, DIF (y e , e~) DAR /^-capture photoabsorption 
a x 10~ 40 cm 2 a x 10~ 42 cm 2 A c x 10 3 s" 1 cr tot x 10~ 26 cm 2 



Closed-shell TDA 

Ohco 3.8 

(0 + l)hu 12.8 

(0 + l + 2)ftw 26.1 

(0 + l + 2 + 3)ftw 31.0 

(0 + l + 2 + 3 + 4)?iw 33.7 

Closed-shell RPA 

Ohuj 1.7 

(0 + l)ftw 6.9 

(0 + 1 + 2)^ 13.2 

(0 + 1 + 2 + 3)hco 16.2 

(0 + l + 2 + 3 + 4)nw 18.2 

Closed-shell RPA + 2p-2h 

Ohco 1.9 

(0 + l)hu 7.1 

(0 + l + 2)hu 13.6 

(0 + l + 2 + 3)ftw 16.1 

(0+1 + 2 + 3 + 4)^ 16.7 

Unrestricted shell model 

Ohio 1.3 

(0 + 7.6 

(0 + l + 2)ftw 11.1 

(0 + 1 + 2 + 3)ftw (13.2) a 

(0 + 1 + 2 + 3 + 4)^ (13.8) a 

CRPA § 17.6(2) 
Expt. 12.4(18) p 



34.6 
39.0 
40.0 
40.7 
40.6 



13.3 
15.8 
20.9 
21.4 
21.9 



14.6 
17.0 
21.1 
21.3 
20.4 



7.1 
10.5 
12.1 

(12.3) a 
(12.5) a 

14.4(1) 



31.0 
62.9 
67.7 
74.0 
75.9 

11.7 
35.1 
41.4 
43.5 
45.4 

13.4 
36.2 
42.5 
43.5 
44.1 



6.4 
31.6 
39.5 
(40.6) a 
(42.2) a 

38.0(7) 



22.7 
29.5 



14.4 
18.3 



17.9 
21.6 



19.6 

(23.6) a 



14.0(12) |2 
14.8(13) pip | 



38.9(9) M 



21(1) U 



Estimate 
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TABLE IV. Inclusive cross sections involving a continuum of 12 N states, but excluding the 
ground state, for various choices of radial wavefunctions: oscillator functions with length parameter 
b = 1.64 fm or 1.82 fm, and Woods-Saxon functions. 

(f^ t , / u _ )DIF (f e ,e~)DAR ^-capture photoabsorption 

a x 10" 40 cm 2 a x 10" 42 cm 2 A c x 10 3 s^ 1 cr to t x 10" 26 cm 2 

Oscillator functions (b = 1.64 fm) 
Closed-shell RPA + 2p-2h 





1.5 


1.1 


1.5 




(o + \)hu 


9.4 


4.9 


22.6 


17.9 


(0 + 1 + 2)hu 


14.2 


5.3 


25.9 




(0 + 1 + 2 + 3)hu 


16.5 


5.1 


24.6 


16.8 


(0 + 1 + 2 + 3 + A)hu 


17.1 


5.5 


25.0 




Jnrestricted shell model 












1.3 


0.3 


0.7 




(0 + \)%UJ 


11.0 


5.6 


27.1 


22.3 


(0 + 1 + 2)huj 


15.9 


5.7 


28.9 




(0 + 1 + 2 + 3)%uj 


(18.5) a 


( 5.5) a 


(27.5) a 


(20.9) a 


(0 + 1 + 2 + 3 + 4)huj 


(19.1) a 


( 5.9) a 


(27.9) a 





Oscillator functions (b = 1.82 fm) 
Closed-shell RPA + 2p-2h 



Ohu 1.6 

(0 + 1)^ 8.9 

(0 + 1 + 2)huj 13.6 

(0 + 1 + 2 + 3)hu 16.3 

(0 + 1 + 2 + 3 + 4)^ 17.2 

Unrestricted shell model 

Ohu 1.3 

(0 + l)hu 10.3 

(0 + 1 + 2)^ 15.1 

(0 + 1 + 2 + 3)hu (18.1) a 

(0 + l + 2 + 3 + 4)/ia; (19.1) a 



1.1 1.6 

5.6 25.5 22.1 
5.9 29.5 

5.7 27.9 20.7 
6.1 28.3 

0.3 0.9 

6.6 30.8 27.5 

6.6 33.0 

( 6.4) a (31.2) a (25.7) a 

( 6.9) a (31.7) a 



Woods-Saxon functions 
Closed-shell RPA + 2p-2h 





0.8 


0.6 


1.3 




(o + i)huj 


6.0 


3.0 


24.1 


17.9 


(0 + 1 + 2)hu 


12.4 


3.1 


27.7 




(0 + 1 + 2 + 3)hu 


14.9 


3.3 


28.6 


21.6 


(0 + 1 + 2 + 3 + A)hu 


15.6 


3.5 


29.9 





Unrestricted shell-model 

Ohco 0.7 0.2 0.8 
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(o + i)hu 


7.0 


3.6 


26.0 


(0 + 1 + 2)hu 


10.5 


3.7 


32.9 


(0 + 1 + 2 + 3)huj 


(12.6) a 


( 3.8) a 


(34.0) a 


(0 + 1 + 2 + 3 + 4)hu 


(13.2) a 


( 4.1) a 


(35.6) a 


CRPA § 


16.9(2) 


5.5(1) 


32.0(7) 






5.1(8) @ 




Expt. 


11.7(18) 


5.7(8) H 


32.8(8) (2. 



19.6 

(23.6) £ 



21(1) E 



a Estimate 

TABLE V. Exclusive cross sections involving the ground state of 12 N only. Woods-Saxon 

radial functions are used, whose asymptotic forms are matched to the experimental ground-state 
separation energies. 

(ffi, DIF (zv e ,e~)DAR ^-capture /3-decay 

a x 10" 40 cm 2 a x 10" 42 cm 2 A c x 10 3 s" 1 B(GT; 1+ -» 0+) 

Closed-shell TDA 

Oftw 2.74 34.5 30.0 2.26 

(0 + l + 2)hLO 2.41 35.4 30.7 2.45 

(0 + 1 + 2 + 3 + 4)huj 2.42 35.2 30.6 2.44 

Closed-shell RPA 

Oku; 1.07 13.2 11.1 0.87 

[0 + l + 2)huj 1.19 18.2 15.1 1.26 

(0 + l + 2 + 3 + 4)?iw 1.24 18.8 15.6 1.30 

Closed-shell RPA + 2p-2h 

Ohuj 1.17 14.0 12.1 0.91 

(0 + 1 + 2)^ 1.16 18.0 14.9 1.23 

(0 + l + 2 + 3 + 4)?iw 1.12 17.0 14.2 1.16 

Unrestricted shell model 

Ohuj 0.59 6.9 5.7 0.45 

(0 + l + 2)hu 0.58 8.4 6.6 0.56 

(0 + l + 2 + 3 + 4)ftw ( 0.56) a ( 7.9) a ( 6.3) a (0.53) a 

CRPA § 0.7(3) 8.9(1) 6.0(3) 



8.9(10) @ 

Expt. 0.66(14) @ 9.1(10) p| 6.1(3) p§ 0.50(3) 



Estimate 
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TABLE VI. Contribution from each multipole to the inclusive cross section in a closed-shell 
RPA + 2p-2h model space. The values on excluding the 12 N ground-state contribution are given 
in brackets. 





(y^n ) DIF 


{v e ,e ) DAR 


/i-capture 




— i n— 40 2 

a x 10 cm 


— i n— 42 2 

a x 10 cm 


A c x 10 J s 


Positive parity 








0+ 


0.11 


0.00 


0.21 


1+ 


2.95 ( 1.79) 


22.52 (1.66) 


15.43 ( 3.87) 


2+ 


2.59 


0.08 


1.36 


3+ 


1.39 


0.03 


0.46 


4+ 


0.66 


0.00 


0.00 


5+ 


0.51 


0.00 


0.00 


6+ 


0.04 


0.00 


0.00 


>7+ 


0.04 


0.00 


0.00 


Sum 


8.29 ( 7.13) 


22.63 (1.77) 


17.48 ( 5.92) 


Negative parity 








(T 


0.07 


0.04 


2.12 


]_- 


3.55 


1.90 


12.25 


2~ 


2.91 


2.36 


7.79 


3" 


1.77 


0.00 


0.11 


4~ 


1.41 


0.00 


0.07 


5 _ 


0.17 


0.00 


0.00 


6~ 


0.17 


0.00 


0.00 


> 7" 


0.00 


0.00 


0.00 


Sum 


10.04 


4.30 


22.34 


Total 


18.33 (17.17) 


26.93 (6.07) 


39.82 (28.26) 
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